Abstract 

In a general class of one dimensional random differential equation the 
convergence of the distribution function of the solution to stationary state 
distribution is studied. In particular it is proved the boundedness respec- 
tively the divergence of the fractional order moments of the solution below 
respectively above some critical exponent. This exponent is computed. In 
particular models it is the heavy tail exponent. When the equation is lin- 
ear this exponent determines a new family of weak topologies (stronger 

(«_^ ' compared to the classical one), related to the convergence to the station- 

Cn , ary state. 
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1 Introduction 

Discrete or continuous time, one dimensional affine stochastic evolution equa- 
tions (ASEE) are studied both in mathematical literature (e.g. Refs.[3[ |T0l[TTl 
[Ml [13]), as well as in physical literature in Refs.fTHl FTTl [T2l [T5] . This interest 
in ASEE comes partly from the occurrence of heavy tail (HT) in the stationary 
probability distribution functions (PDF), in the models of physical, economical 
and biological processes. The ASEE are also related to reduced models of the 
self-organized criticality phenomenology in plasma physics fRef.|15j). and to the 
renewal processes (Refs.[3lfT0l [11]). 

The connection of the our approach to the linear ASEE models with HT 
in the stationary PDF is the following. Denote by X-t:{uj) a solution of ASEE 
and pt{x) :— probi^{\Xt{uj)\ > x). Suppose also that the ASEE has a stationary 
PDF, so we denote p{x) = \im pt{x). We define the heavy tail exponent ;9c 

by the asymptotic estimate p{x) — 0{x^^''l{x)), where l{x) is a slowly varying 
function. 

The occurrence of the HT in the stationary PDF is related to the dynamical 
effect (Refs.[13j [15]) that will be studied in this article. When the stationary 
PDF of the solution Xt{u}) of the ASEE has HT with exponent /3c, then for 
i — >■ oo the fractional order moments EflXjl^] remains bounded for < p < /3c, 



respectively diverges for p > (3c- This is related to the "variance explosion" 
phenomena studied recently in the mathematical finance (Ref.[13]). In this 
article we obtain a simple, explicit formula for (3c- 

Previously explicit algebraic methods for computing (3c were elaborated in 
some special cases: in the framework of the discrete time models in Refs. [21 HZ], 
with i.i.d. additive and multiplicative noise, and in the case of i.i.d. additive 
noise and multiplicative noise modelled by a finite state Markov process in 

Ref.[n]- 

In the continuos time case, with the multiplicative noise modelled by a finite 
state Markov process, rigorous foundation of the computation oi (3c was obtained 
in Ref.[10_. In Ref.[T5] the multiplicative and additive noise were modelled by 
a superposition of Ornstein-Uhlenbeck processes. An explicit formula for (3c 
was derived by asymptotic methods. In all of the cases the exponent (3c is 
independent of the additive term. 

The ASEE model equation that is considered here is a class of one dimen- 
sional random differential equation (RDE ), which extends previous results from 
Ref . [T5] ■ by using new topological vector space methods. The additive and the 
multiplicative random terms in our model are stationary processes. The multi- 
plicative term is a generalization of the stationary Gaussian processes, having a 
very general, possibly algebraic, correlation decay. Our results and those from 
Ref. |10) are complementary. 

The main part of the new results are related to the convergence to stationary 
state in a class of linear RDE. Nevertheless, the formula for (3c was extended to 
a class of non linear RDE. 

Define the subspace of real valued continuos functions Cy{M.) C C(R) by 
the condition / S Cj{R) <;=> |/(a;)| — o{x~'^)- For the class of linear models 
considered here, we prove the existence of the suitable defined weak limit of the 
PDF of Xt, in the sense that there exists a random variable Xoo, independent of 
the initial conditions, such that if Xt is the solution of the RDE and < 7 < /3c, 
then limE[/(Xt)] = E[/(Aoo)] for aU / e C^(R). The topology on the space 

t—>oo 

of probability measures, related to this new class of weak limit, is stronger 
compared to the classical weak topology. 

In particular we have f{x) = \x + z\p e C7(M) for all z G C when < p < 
7 < (3c, so for i ^ 00, we have E [\Xt + z\p] -^ E [\Xoo + z\p]- Moreover we 
prove that if p > /3c then E [\Xt + z\p] -^ 00, for all but eventually one initial 
conditions. The expectation values, E [\Xt + z\p] , E [|Xoo + z\p] are related to 
a class of generalized LP spaces that includes also the p G ]0, 1[ case. The proofs 
uses the geometry and topology of these general L^ spaces. In particular we 
prove that X^o G L^ for all < p < (3c- When heavy tail exists this property 
determines the heavy tail exponent (3c- 

We apply these L^ methods to the study of nonlinear RDE , when the non- 
linear term has a weak nonlinearity. In this nonlinear model we derive a weaker 
results. There exists also a critical exponent (3c, such that for < p < /3c the 
moments E[|Xt -I-2:|p] are bounded for large t , for all initial conditions. For 
p > (3c and for sufficiently large initial conditions the moments E[|X( + z\p] 



diverges exponentially for large time. The same formula for /3c give the demar- 
cation line between the order of the divergent, respectively bounded fractional 
moments. 

In all of the cases /3c is given by a simple analytic formula in term of 
two physically significant parameters of the multiplicative term only, suggesting 
some non commutative Central Limit Theorem on the afRne group. 

The structure of this article is the following. Section [2] provides the nota- 
tions, the linear model and the main theorem. Section [3] gives a sequence of 
propositions leading to the proof of the main results. The study of the non lin- 
ear model, by using results from Section[3]is presented in Sectional A summary 
of the results is given in Section [SJ 

2 Description of the linear model and the re- 
sults. 

2.1 Notations and definitions. 

The stochastic processes are defined in a fixed probability space {O, T , P} 
with expectation value E(^[/(a;)] = E[/] = J^f{uj)dP{uj). By uj will be de- 
noted a generic element of ft. Two driving J^— measurable stochastic processes 
{Q{uj), 4>t{i^) : K X fi — >■ R^, and a constant a > (the instability threshold) 
defines formally the linear RDE : 

^5^ = -(a + C*M)x,H-t-0,H (1) 

whose solution can be constructed explicitly (see below). The notations 
Xt{uj) or Xt, for the solutions of RDE will be reserved. Without loss of generality 
we impose E[^t] = and deterministic initial conditions Xo(w) = xq. The 
argument uj will be omitted when no confusion arises. We shall denote by S the 
set of solutions of Equation ([Ij . 

We will preserve notation L^ also for some unusual Lebesgue spaces with < 
p < 1 (see Ref.0, page 75): 

Notation 1 Let p > and denote ap ;= min(l,p). We define ||/(w)|L := 

(E[|/|P])''''/^ and LP := LP{Vl,T,P) = {/|||/||p < oo}. More explicitly for 

p > 1 we have the usual norm ||/(w)|L := (E[|/|P]) , while for < p < 1 we 
have the distance to the origin given by ||/(i^)|L := E[|/|'']. 

These complete metric vector spaces were already used in the study of prob- 
ability distributions having HT (Ref.[7]). 

Remark 2 Also for p G (0, 1) the topology induced by the distance \\f — g\\ 
and the vector space structures are compatible. This follows from the general 
inequality 

\\afiu)+giu)l<\a\'^'^ ll/llp + llffll, (2) 



where p > and f{uj), g{uj) G L^. All of the LP spaces are complete (see page 
15 from Ref.J^). It is easy to see that for p € (0,1) the unit ball is not convex 
and in nonatomic cases their dual is trivial. 



When < p < 1 the Inequahty © resuhs from \a + bf < \af + \bf . 

The final convergence resuhs will be formulated in the terms of a subspace 
C^(M) of continuos functions on M. 

Definition 3 For any 7 > we define the subspace C-y(R) of the space of the 
continuous functions C(M) by the condition: 



l/WI 

|a;tToo(l + \x\y 



f e C^(R) ^ , hm ^/\j' = 



and the topology on C-),(R) by the norm 

In the proofs the extension of the space of Holder-continuous functions to 
the whole real line proves useful, in the study of the tail effects with extreme 
delocalization, i.e. /3c ^ 1: 

Definition 4 If a e]0, 1] then f{x) e Ha ^ {=l(ci > 0) such that \/{x,y G 
K) [\fix + y)-fix)\<c,\y\'-]}. 



We will preserve the same notation, to define a class of functions useful in 
the study of the more localized regimes, when the stationary PDF has /3c > 1: 

Definition 5 If a > 1 then f{x) e Ha ^ {3{gi{x) £ Hi, q e C ) s.t. f{x) = 



> 



Remark 6 If f{x) G Ha then \f{x)\ < a + 6|x|", both for a = I, f( 



< 



or some 



constants a,b. The exact values of the constants Ci,a,b are irrelevant. 



Remark 7 It is elementary to check that if z (z <C , < a, then f{x) 

\z + xT G Ha- 



By using the spaces Ha, C-y(M) we define a parametrized families of weak 
topologies on the set of probabilistic Borel measures on M. These topologies, 
with indexed by a fixed p > 0, are defined as follows: 



Definition 8 Denote JCp := ^ Tia (finite linear combinations) . For t — >■ 

Qe]o,p[ 
oo, the sequence /it of the Borelian measures on M converges to the stationary 
measure jjLoo in the topology T{ICp) iff V/(x) G /Cp we have lim L f{x) d^,t{x) = 

J^f{x)dfioo{x) 

In the final formulation of the results, the asymmetry in the definition of 
the spaces T-ia will be removed. An apparently stronger class of topologies are 
given by the following 

Definition 9 For t -^ cxd, the sequence /if of the Borelian measures on M con- 
verges to the stationary measure fio^ in the topology T{C^) iff V/(x) G C^(R) 
we have lim ^^ f (x) dpLt{x) — ^^ f {x) dpLoo{x) . 

We observe that the T{G^ ) is strictly stronger than the weak topology. 

2.2 Specification of the linear model. 

Without loss of generality and for sake of simplicity, we suppose that a suitable 
rescaling of the time variable was performed. Consequently, the multiplicative 
noise C,t obeys the following 

Condition 10 The stochastic process Ct is centered and stationary. There ex- 
ists a unique "diffusion constant" D > and a set of positive constants K'r , all 
independent of the times t and s, such that for p > Q, s > \ and s — 1 <t < s 
we have: 

Kp < exp {^Dp^ s) E [cxp {(1 -p)Ys- Yt}] < K+ (4) 

where 

Yt := / C. ds (5) 

Jo 

The exact values of the constants K^ and K^ are irrelevant. This Condition 
is satisfied by a large class of Gaussian processes, as shown in Proposition [Ml 

We will prove that the large time asymptotic behavior of the solution of the 
RDE (H)) is determined by the parameter a from Equation ([Ij and the constant 
D from Equation (U). So we introduce the following 

Notation 11 ; The main results are expressed in terms of the critical exponent 

/3c := a/D (6) 

Another important quantity is ^p :— Gp D {p~/3c) withp > 0, whose sign changes 
atp^ Pc- 

On the additive noise (j)t{'^) we impose the following 

Condition 12 1. The process 4't{'^) is stationary. 



2. The processes Ct{^), 4't{^) o,fe independent. 

3. We have ¥.[\(j)t{uj)f] < rUp < oo, \/p e]0, /3i[, where (3i > niax(l,/3c). 

4-. Reversibility: the equality in distribution (pti^) = 't>-t{^)- 

From Conditions [T2l item |3] results that the additive noise is allowed to have 
a HT with exponent larger than max(l,/3c)- 



2.3 The results, for the hnear model. 

Under the previous notations, Conditions [TU] and [T^J we have the following 
results, that will be demonstrated in the next sections. First we have the main 
Theorem and its Corollary 

Theorem 13 Let } e C^iji.) andO < j < f3c. There exists Xoo{uj) G flpelo /3c[^^ 
such that: 

1. fiXoo) e L^ 

2. We have \imE[f{Xt)] = E[f{Xoo)]- In particular if p e]0,/3c[ and z G C, 

then lim E \Xt + z|^ = E \Xoo + zf < oo. 

3. If p > /3c then, except for at most a special value of the initial condition 
Xt=a, we have lim E \Xt + z\^ = oo. 

This Theorem leads immediately to the following Corollary: 

Corollary 14 Denoting by Ft{x) and Foo{x) the cumulative PDF of Xt and of 

Xoo, respectively, then: 

1. If-f e]0,/3c[ and f e C^iR), then \imJ^f{x)dFtix) = J^f{x)dFooix). 
In particular, for any complex constant z we have lim Jj^ |x + 2:|^ dFt{x) = 
/jj \x + z\'^ dFoo{x) < oo. 

2. If ^ > Pc then except for at most a special value of the initial condition, 
we have lim f„ \x -\- z\'' dFtix) = +oo. 

The Item 1 of this Corollary states that weak convergence of the measure 
generated by Ft{x) in the topology T{Cj). 

The previous theorem results from the following technical Lemma and its 
Corollary 

Lemma 15 There exists Xoo{^) G Hpelo B \^^ such that: 
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2. { fix) 6 ICpJ ^ { \imE[f{Xt)] = E[/(Xoo)]}. In particular forp e]0,/3. 



and z e C, we have lim E \Xt + zf =E \X^ 



3. li p > Pc then, except for at most a special value of the initial condition 



^t=o, we have lim E \Xt + z\ — oo. 



This Lemma leads immediately to the following Corollary: 

Corollary 16 Denoting by Ft{x) and Foo{x) the cumulative PDF of Xt and of 

Xao, respectively, then: 

1. If p e]0,^c[ and f(x) e Up, then lim J^f{x)dFt{x) = J^f{x)dFoo{x). 
In particular for any complex constant z we have lim L \x + z\f dFt{x) = 
/g \x + z\P dFocix) < oo. 

2. If p > (3c then except for at most an eventual special value of the initial 
condition we have lim f„ la; + z\p dFt(x) = +oo. 

The Item 1 of this Corollary states that weak convergence of the measure 
generated by Ft{x) in the topology T{JCp). 

We will prove that the topologies T{JCp) are equivalent T{C^), when re- 
stricted to probabilistic measures. 

Remark 17 The convergence/ divergence of the moments for p ^ /3c does not 
imply the existence of the HT. Indeed, if(f)t(j-o) — andp e ]0, f3c[, then lim ||^t|L 

\\Xoo\\p = 0. See Remark\M 

Remark 18 From the proofs will be clear that for p < (3c, when 7^ < (see 
Notation \ll\) the speed of the convergence for large time o/||Xt|| is 0(exp(7j,i)). 
If p < (3c then the speed of divergence is ||^t|L — 0(exp(7pi)). 



3 Proof of the Theorem [13 



For technical reasons we introduce the following additional notations 
Notation 19 

At :— exp {—at — Yt) 

Bt := / (t)rAt/Ar dr 
Jo 

Ht := I (prAr dr 
Jo 

d{t)^E[Yt']/2 

Furthermore, all the constants denoted by Kn with n integer, are independent 
of the time variables t, s, x, y. 



3.1 Some initial results. 

Because the driving noises are classical functions, then according to the Nota- 
tions [12] the rigorous, integral form of Equation ([l}, with the initial condition 
Xo, is 

Xt ^xoAt + Bt. (7) 

Remark 20 From Equation Qj and Equation 0) we get: 



\xor'\\Ai~\\Bt\\^ 



^wxti^ixor^wAii + WBtW^- (8) 



The following proposition is the key point in the difficult part of the proofs 
that needs L^ bounds when < p < 1. This case include also the study of the 
"extreme heavy tail" effects (Ref.|15)V when /3c is very small. 

Proposition 21 Let U{uj) > and V{ijj) > be random variables such that 
V, UV e L'^andE[V]>0. Then for p £]0,1] we have 

E [UP V] < (E [U V]f (E [V])^~P . (9) 

Proof. This result follows from Jensen inequality applied to the concave func- 
tion X — > xP. Define a new expectation value Ei [/] :— E [fV] /E [V]. 

For any concave function g{x) : R — >M by Jensen inequality we have Ei [g{f)] < 
(/{El [/]}. Consider now g(x) :— x^ and /(w) — U{uj). We obtain 

E [UP V] /E [V] < (E [U V] /E [V] )p 

that leads to the Inequality (0). ■ 

The following Proposition results from Condition [TOl bv setting s = t: 

Proposition 22 Under the Condition ] 1 (A there exist positive constants D, K^ , 
K~ , Ki, and K2 such that \fp G ]0, oo[ and t > we have 

K- exp(p2 Dt)<¥. [exp(-p Yf )] < K+ cxp(p2 D t) (10) 

K2 exp(t7p) < Ptilp < ifiexp(t7p) (11) 

Remark 23 From Inequality \11\) and Notation ] 11\ results that for < p < /3c 
we have ||A(|| — >■ when t ~¥ 00. In the particular case (j)t{^) — 0, by Equation 
^ we obtain also [[Xt[[ -^ 0. 

Proposition 24 Suppose the process Q is Gaussian, stationary, centered, with 
continuous realizations and has the correlation decay for t ^- 00 

C{\t[) = E{Ct+r Cr) = 0{1 + [tr'-'). (12) 

Then Condition \10\ is satisfied and the Inequalities ^ \1(A 111]) results. 



Proof. Because (t is stationary, we will use a method similar to the derivation of 
the Taylor-Green-Kubo formula (Refs.[2l[T8l|4l[5]). For the sake of completeness 
we adapt the proof By using Equations ©[H]) and the Notations [121 we obtain 

t 

2 d{t) = E[Yt^] =2 f{t- x) Cix) dx (13) 



Combined with Equation P^ it results that for t — ;> oo we have d{t) = Dt + 
0(1), where 

/•OO /"OO 

D= C{\T\)dT= E[C.Co]dT (14) 

Jo JQ 

Also from Equations ([T2l[T3)) for |a;| < 1 we have 

d{t + x) - d{t) =^ 0{l) (15) 

As a consequence, there exists a constant K^ such that for all f > we have 
uniformly in t: {Dt- K3) < d{t) := ¥.[Y^]/2 <{Dt + Kz). 

Because Yt has stationary increments we obtain VJ\Yy Yt] = d{t) + d{y) — 
d{y - 1). With Z = {I - p) Ys~Yt we get E [Z'^/2] = p^ d{s) + [p {d{t) - d{s)) + 
{l-p)d{s-t)]. 

From \s — t\ < 1 and Equation (|15p it results that the last bracket is uniformly 
bounded. Since Z is Gaussian, we obtain the Inequality Q and by Proposition 
[221 the Inequalities (fTOl lTT I) . ■ 



Observe that by Equation p4l) the constant D is related to the zero frequency 
component of the correlation function. We shall prove the following general 

Proposition 25 Let p e]0,1[, < y — I < x < y and suppose that 

1. The stochastic processes Yt,Ct satisfies Condition \l(A 

2. The stochastic process ft is stationary, independent of Yt and we have the 
bound uniformly in t: 

|E[/t]|<if4<oo. (16) 



Denote 

rv 



ft At dt 



(17) 



bf{x,y) := 

Then for some constant K^ we have 

bfix,y) < K^expiyjp) . (18) 

Proof. Since for p e]0,l[ we have ap — p and \\g\\ = E,[\g\'^], the Equation 
([TTI) can be rewritten as bf{x, y) =E [U'' V], where 

V{uj)^exp[-pYy{L,)] 

U{uj) = / ft{uj)At{uj) exp[yj,(w)] dt. 
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Recalling Proposition ?IT\ and Inequality (O we get: 

bf{x,y)<{E[UV]rmv])''^ . (19) 



Then, on the one hand, according to Proposition [221 and Inequality (|10l) . the 
term E [V] in Equation (fT9|) is bounded by 

E[V]<Keexp{p^yD) . (20) 

On the other hand, due to the independence of ft and Yt we obtain 

E[UV]= I E[f]E[Atey.^{{l-p)Yy}]dt. 

Recalling Inequality (ITCl) , and the definition of At , it follows that 

\E[UV]\<Ki I exp(-ai)E[exp{(l -p)y2,-rt}] dt 



which, with the help of Inequality ([4]), reduces to 

\E[UV]\<K'jeyi^{-ay+p^Dy). (21) 

After simple calculations, recalling Inequalities (IT^ EOl [^T|) . the Inequality (|18p 
is obtained. ■ 

Remark 26 From the Condition llSl for fixed T, by symmetry and stationarity, 
we have an equality in distribution of the direct and reflected process (pt = (f>T-t- 

From the previous Remark, it follows the following 

Proposition 27 // f{x) E Hp then K[f{Bt)] = K[f{Ht)] and in particular 
\\Bt + z||p = \\Ht + z||p, where z e C, p > 0. 

Proof. From the Definitions |1][5] results, that if f{x) 6 Up then f{x) is contin- 
uos, hence measurable. From Remark[5]we have |/(a;)| < C1+C2 [a::]''. According 
to the Notations [T51 

Bt^ ( (t)r exp {-a{t - r) - {Yt - Yr)) dr 
Jo 

From Condition [TUl results that the process Y^. has stationary increments, so for 
fixed t we have the equality in distribution of the processes Yt —Y^ = Yt-r- On 
the other hand, from Condition [5] results that because the processes (PtiYt are 
independent, we have the equality in distribution 

{^r,Yt-Yr) = {4>r.Yt-r) 
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Consequently 



Bt= 4>r exp (-a(t - r) - Yt-r) dr 

JQ 



By the change of the integration variable t -^ t — t and from Remark [ 
results the equality in distribution 



Bt = (pt-T exp {-ar - Yr) dr 
Jo 



(22) 



According to the Remark I26[ for fixed t we obtain (/>,- ~ (j)t~T- By the inde- 
pendence of the processes 0r,Crj consequently the independence of (l)T,Yr we 
have 

{</),, Yr) ^ i(bt-r,Yr) 

So from Equation (P^ results 

/■* 
Bt — 4>T exp (— ar — 1^) dr 

or (see Notation [TO]) St = Ht which completes the proof. ■ 

The following Lemma and its Corollarv 1301 will be used in subsequent works. 

Lemma 28 Under the previous Notations \11[ I jffi Conditions [73 and \12[ we 

have 

1- IIv>u>0,^p^Q and < p < /3i then we have the L^ hound 

\\Hy - i/„||p < i^8 {exp(u7p) +exp(«7p)} (23) 

for some constant K^ independent of u,v. 



2. There exists Xqo G no<o<;g^ ^'' such that Jim \\Hf — Xoc||p = when 

pe]0,/3J. 

3. Vp e]0,/3i[ we have 



pe]0,/3J. 



||St||p<i^8[l+exp(i7p)] (24) 

Proof. 

1. If pe [1, /3i[ t/ien (see Notation[Tll) we obtain ||iJ„ - Hu\\p < /J ||(?!)s A^Hp ds. 
Recalling Conditions [T2l items 2 and 3, it follows that 

Consequently, by using the Inequality pT]) results 



\Hv - Hu\\p< rripKi / exp {t 7p) ds 
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which immediately lead to Equation 

In the "very heavy tail" case when p g]0, 1[ we use the Proposition!^ 
where we replace ft = (t>t and Equation p7)) . We introduce the notations 
[x] for the integer part of x, the notations n^ := [w — u], w_ :— u + n^ and 
use the decomposition 



brArdr =(/ + h/ +/ ) (l)rArdT 



We get 



fe=i 

By using Proposition [511 Equation ([TS]) with ft = (j)t, after simple estima- 
tions we obtain Equation (j23p . which completes the proof. 

2. Let p e]0,/3c[. Then we have 7p < and from Equation (P5|) it results 
that if t„ ^^ cxD then i/t^ is a Cauchy sequence in L^. 

The LP spaces, including p e]0, 1[, are complete (Ref-IS], page 75), there- 
fore an Xoo(w) G LP exists such that \\Ht — X^oWp — > 0. Because < p < 
q ^ LP Z) L'', it results that Xooioj) £ r\Q<q<p, -^'• 

3. It results from Proposition[27land Inequality (l23l) . by setting u = 0, v = t. 



Lemma 29 Let f{x) e T^p /ixed. Then "^ e LP ^ f{^) e L^ .The correspond- 
ing application LP — > C, defined as i^ 9 ^ — )• E [/(^)] G C, is continuous. 

Proof. From f{x) e T-Lp results |/(a;)| < a + 6|x|P (see Remark [5]). Hence 
^ G LP ^ /(*) G L^. 

For the proof of the continuity we consider first the case < p < 1. Let 
^,X^ LP and recall that in this case f{x) £Hp ^ \f{x + y) — f{x)\ < ci \y\P. 
Then it follows that 

IE [/(* + x)] - E [/W]| < E !/(* + x) - /(*)| < ciE Ixr = ci llxllp 

which proves the continuity in the case of the (unusual) LP space. 

In the case p > 1 we observe that according to the Definition [S] it is sufficient 
to prove the continuity for the set of functions of the form f{x) = \g{x)\^ 
when g{x) G Tii, that generates the space Up by finite linear combinations. 

Let ^, X G LP with p > I and recall that in this case f{x) G TLp ^ 
\g{x + y) — g{x)\ < ci \y\. So it is sufficient to prove that 

||x||p^O^E[/(vI/ + x)]^E[/(^P)] (25) 
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In the case f{x) = \g{x)\^ the previous Equation (l?5]) is equivalent to the 
condition: 

11x11,^0 ^||,9(* + x)ll,^ll3WII, (26) 

Observe that from the usual {p > 1) Holder inequality and from g{x) 6 Hi we 
obtain 

M^ + x)l Mm,\ < Il3(* + x)- .g(*)ll, < llcixll, 

which completes the proof for the case p > I. ■ 

Corollary 30 Let p e]0,^c[ and f{x) e Hp. Then IE[/(i7t)] -> E[f{Xoo) and 
E[f{Bt)] ^ E[f{Xoo)] for t -^ oo. 

Proof. The convergence of E[/(i?t) results from Lemma [^ part [^ and Lemma 
m From Proposition [m resuhs E[/(Bt)] = E[/(iJ()], which completes the 
proof. ■ 

Remark 31 The distribution of X^o G L^ , luith p £]0,/3c[ is identical with 
limt^oo Ht :— L (prArdr, where the limit is in L^ . Clearly, in generic cases 
Xao is non degenerate, and for the limiting measure i-i{x) — prob{\XQo{^)\ > x), 
we obtain < J^ y^d^ (y) < oo for all p g]0, /3c[- 

3.2 Proof of the Lemma I15L 

Proof. The random variable Xca G no<p<^^i^ was constructed in Lemma [^ 
part[2l Let p g]0,/3c[ and f{x) G Hp. If f{x) S /C^^, then it can be represented 
in the form f{x) — Y^'^^i fi{x) where fi{x) € Hp- for pi e ]0, /3c[ and I < i < n. 
It is therefore sufficient to consider the case when f{x) G Hp for p g]0,/3c[- 

1. It is sufficient to prove that {fix) € Hp,p e]0,/3c[} => {fiXoa) G L'^}- 
We now use X^o G L^ and Lemma [^ 

2. Let t — >■ oo. It is sufficient to prove that {f{x) e Hp,p g]0,/3c[} ^ 
{ limE[/(Xt)] = E[f{Xoo)]}- From Proposition [27] it results that 

E[fiBt)]=nfiHt)]. 
Consequently, according to Corollary [30l we have 

limE[/(i3,)]-E[/(Xoo)]<oo. 

We will use now Lemma [29] with x — xo At and '^ — Bt- Because for < 
p < j3c^e. have exp(t7p) -^ then, by Inequahty (fTTj) results ||^t|| — > 0, 
so llxIL ~^ 0- Consequently by Lemma [29] we obtain 

E[f{Xt)] = E[f{Bt + xa At)] ^ E[f{X^). 

E[f{Xt)] = E[f{Bt + xoAt)] -^ E[f(Xoo)- The last remark follows from 
the fact that f(x) G Hp for all z G C if the function f{x) is defined by 

f{x):=\x + z\P. 
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3. Consider the case when p > (3c i.e. jp > and therefore exp(7pi) -^ +oo. 
Suppose, ad absurdum, that there exist two initial conditions xq and x'q 
such that the corresponding solutions Xt{uj) and Xl.{uj) of Equation ([T]) 
are bounded in L^. 

The stochastic process Zt{uj) — Xt{uj) — X'f.{uj) is also bounded and it is 
determined by the homogenous equation. According to Equation ([7]), it 
results that 

Ztiio) = (a;o -Xq) Af(tj) 

and thus \\Zt\\ — \{xo — x'q)]'^" ||^t|L- Then, recalling the first part of 
Inequality (|lll) we finally obtain 

\\Zt\l>\{xo-x'X' K2exp{^pt) (27) 

Because when p > f3c we have, jp > and therefore exp(7pt) — >■ +oo., by 
Equation ^7} the proof is completed. 



3.3 Proof of the Theorem [HI 

The main point of the proof is the density of the space Hq in the space C^ (R) . 
For the convenience of the reader, we recall some very general definitions and 
a generalization by Nachbin of the Kakutani-Stone density theorem from (Ref. 
[8], Theorem 2 ) to the non compact case. See also (Ref. [6]). 

3.3.1 The Nachbin approximation Lemma 

We use the terminology and notations from (Ref. [6]). Let X a completely 
regular topological space. By C{X) we denote the space of continuous real 
valued functions on X. 

A directed set V, where V C C{X), is a set such that for any i'i,W2 S V 
there exist A > and v &V such that vi{x) < Xv{x) and V2{x) < Xv{x). 

We suppose also that for any x G X there exists v E V such that v{x) > 
(i.e. the set V is pointwise strictly positive). A function /(.) G C{X) is an 
element of the weighted space CVoo(X) if and only if for all v G V the function 
v{x)f{x) vanishes at infinity. 

The topology on CVoo{X) is defined by the seminorms indexed by elements 

of V, denoted Pv{f) 

Pvif) = sup \v{x)f{x)\ 
xex 

where f eCVooiX), veV. 
Recall, a lattice L C C{X) is a set of functions closed under min and max 

f,g<EL^ min [f{x),g{x)] G L and max [f {x) , g{x)] G L 

The following Lemma of Nachbin (Refs [SlIS]) will be used 
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Lemma 32 (Nachbin JS^) Let X he a completely regular space, V a pointwise 
strictly positive set of weights, L a sublattice of CVoo{X) and f € CVoo{X). 
Then f can be approximated by elements of L in the C'VooiX) topology if and 
only if for any x,y E X and e > there exists g E L such that 

\f{x)-g{x)\<e (28) 

\f{y)-9iy)\<£ (29) 

The following Corollary results 

Corollary 33 For any < a < 7 the space TLa is dense in C^(R) in the 
topology induced by the norm from Equation (0). 

Proof. In order to apply the Lemma |32] in the our case, we consider X = R and 
the set V consists of a single function v{x) = (1 + la:])"'''. In this case the space 
CVoo{X) is C-y(R) from Definition [31 First we recall RemarklHland observe that 
'Ha is contained in C-y(R) when < a < 7. 

To apply Lemma [321 for any fixed a > we define the subset L of 7ia as 
follows: h € L ii and only if there exists 77 G ^^l such that for any x G R we 
have h{x) = |77(x)|". It is easy to see that in this case h S Ha for any a > 0. 

Indeed, for a > 1 it follows from the very definition of Ha- For < a < 1 
we have ||?7(a;)| — \il{y)\ I ^ \v{^) ~ viy)\ ^^d now we use that rj G Hi, i.e. 
\t]{x) - ri{y)\ < c\x~y\. So we obtain \h(x) - h{y)\ = \\v{x)\°' - \viy)\°'\ < 
[c\x-y\]" 

It is easy to verify that the set L defined below is a lattice. 

Indeed, suppose that /ii,/i2 G L C Ha, for some fixed a > 0. This 
means that they can be represented as hi{x) — |77i(x)|", h2{x) = |772(a;)|", 
with r]i^2 £ Hi- Then we have h{x) — min [/ii(x), ft.2(a^)] ~ 1^7(2^)1" j where 
r]{x) = min[r]i{x),r]2{x)] £ Hi, because Hi is closed under min, max opera- 
tions. In a similar manner we prove that L is closed under max[]. 

In order to use the Lemma [32l it remains to prove the restrictions given by 
Equations (|28l[29l) . 

Moreover, because the positive and negative part of a function from C.y{M.) 
belongs to C^(R) too, it is sufficient to verify that any non negative function 
f{x) > from C^(R) can be approximated by elements of L. 

But it is easy to check that by using the family of functions g{x) = \ax + 5|" G 
L, a,b €C , the conditions given by Equations (|^ [^ of the previous Lemma 
[221 are obeyed, for f{x) > 0, which completes the proof. ■ 

3.3.2 Proof of the convergence of the measures (Theorem I13p 

Proof. We will use Lemma [151 and its Corollary [TBI In order to complete the 
proof, by extending the results from Lemma [TSl and Corollary \T^ to Theorem 
[T51 it is sufficient to prove that form the convergence of the measures in the 
weak topology T{lCp^ ) it follows the convergence of the measures in all of the 
topologies T{C^) where 7 < /?c . 
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In this end it is sufficient to prove that from hm /^ /(x) dFt (x) — J^ f{x) dFoc {x) 

for f{x) £ Ha, a. g]0,/3c[, where Ft[x),Faa{x) where defined in Corollary[Tni it 
follows that for any g(.) G C'7(K), < a < 7 < /3c we have the same convergence 

lim / g{x)dFt{x)= / g{x) dFoo{x) 

*^°°JR JR 

Denote ut {x) :— Ft{x) — Foo{x) . From Corollarv [TBI results that there exists 
some finite constant K^ such that 



{1 + \x\y \dutix)\ =X^, 0<7</3c 



(30) 



We emphasize that this is the point where ^ < l3c is used. Let g{-) & C^(R), 
and according to the previous Corollary [33] consider the sequence /„ (x) G Ha 
of approximants g{x) in the topology of C^(R) where < a < 7 < /3c- So we 
have 

9{x) - fn{x) 



lim sup 



ii + \x\r 



limp^(g-/„) =0 



(31) 



We have to prove that for any e > there exists T^ such that it t > T^ then 
\Jg^g{x)di/t {x)\ < £■ We have 

g{x)d,.t {x) = f U{x)dut (x) + f g(^)-/"y (1 + 1^1)7 d^^ (^) 

Jr Jr (1 + fI) 



'j{x)dvt {x) 



< 



fn{x)dvt {x) 



R 



+ P7(5-/n) / {l+\x\y\dvt{x)\ 



R 



By Equations ([301 [3T|) we select n such that p~^{g -~ fn) < ^/i'^K^). 

Because fnix) G Ha, a s]0,/3c[ from Corollarv [TBI results that there exists 
Te such that if t > T^ then \l^ fn{x)di^t ix)\ < e/2 . 

This last inequality completes the proof of the convergence in the topology 
T{Cj) , ^ < j3c, under the hypothesis that we have convergence in the topology 

Consequently the results from Lemma [15] and its Corollary [TH] can be ex- 
tended to Theorem [13] and Corollary [14] which completes the proof. ■ 

4 Application 

We consider now a generalization of the Equation ([T|), containing a new non 
linear term ^t [Xt{uj)^uj]. We write the equation symbolically as 



dXtJL^) 
dt 



-(a + Ct(c^))^t(^) + *t[^t(w),w] 



(32) 



The driving noise terms Ct{'^), ^t [Xt{'^), i^] satisfy the following restrictions. 
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Condition 34 The stochastic process Ct('^) obeys Condition \l(A There exists a 
stochastic process (/>t(w) > such that 



\^t[Xtiuj),uj]\ <(t)t{u) a.e. 
and (j)t{io) obeys Condition M'A 



(33) 



In this case the existence and uniqueness of the solution is guaranteed by 
Theorem 221, page 65 from (Ref. [T]). The heavy tail effect is manifest now by 
the boundedness/divergence of the moments E [|Xt + z|''] when p ^ (3^- 

Theorem 35 Denote by Xt{uj) the solution of the Equation i32\) . Under Con- 
dition^^ we have the following behavior for t ^ oo. 

1. when < p < (3c the fractional moment E[|Xi + z|''] are uniformly 
bounded in t. 

2 . For p > Pc and for sufficiently large values of Xq , the moments E [ | Xf + z\^] 
diverges. 

Proof. In analogy to the Equation ([7]) and Notations [T9l we have the rigorous, 
implicit, integral form of the Equation ([5^ 



Xt =xoAt + Bt 

%:= [ 'i't[Xtiuj),uj]At/ArdT 
Jo 



(34) 
(35) 



Because ||/(a;)||p :== {E[\f\P]f''^P (see Notation [T]) it is sufficient to study 
the boundedness of \\Xt + z\\ , and by Equation ([2]), the boundedness of \\Xt\\ . 
From Equations (p4l [2]) results 



i^or^i^di. 



Bt 



<II^*Ilp<|-T^0| 



\At\ 



Bt 



According to the Inequalities ([TTl [551) , we have 



\xq\'''' K2 exp(i7p) 



Bt 



< \\Xt\\p< Ixol"" Kiexp{t-fp) + 



Bt 



(36) 



(37) 



So it is sufficient to prove (in analogy to Inequality 



the new bound 



Bt 



<Ks [l + exp(i7p)] 



Observe that because At/A^ > 0, from Equations 
inequality 



Bt{uj) 



< I 4>t{^)At/ArdT 




(38) 
we have the 

(39) 
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According to the Condition [Ml in the Equation ([5^ we can identify 

/ (t)t[uj)At/ArdT^Bt 

Jo 



.FinaUy, by using Inequahty (|24)l we obtain the Inequahty ([38)) . From Inequal- 
ities d^EHl), it resuhs 



respectively 



XtWp < kor'' Ki exp (t7p) + Ks[l + exp{tjp)] (40) 



\x0rK2 exp{tjp)-Ks [1+eMtlp)] < \\Xt\L (41) 



If p < /3c then we have 7p < and from Inequality (|40| results that || Xt|| is 
bounded. If p > /3c and if the initial conditions are sufficiently large, such that 
\xo\'^'' > Kg/K2, then we have 7^ > so by Inequality (|1T|) results that || -^^tL 
diverges exponentially, which completes the proof. ■ 

5 Conclusions. 

In a class of one dimensional random differential equations the large time behav- 
ior of the solution was studied. When the equation is linear we proved that the 
convergence to stationary state can be described in the framework of new class 
of weak topologies on the set of probability distribution of the solution. These 
topologies are stronger, compared to the classical weak topology. The study of 
the weak convergence involves in a natural way the study of the convergence of 
the fractional order moments of the solution. A critical exponent /3c was defined 
such that the moments of the solution, of order p remains bounded if p < /3c and 
diverges, on a massive set of initial conditions, when p < f3c- When heavy tail 
exists then /3c is the heavy tail exponent. The speed of convergence/divergence, 
for large time, of the moments of order p of the solution is exponential (see 
Remark [18)) , depending on p — (ic 

The strength of the new family of topologies, that describe the approach to 
the steady state, increases with /3c. 

An important result is the exact and simple Equation ^ for /3c, in term 
of the parameters from the multiplicative term only. The convergence in this 
topology of the distribution function to the stationary distribution was proved. 
A new topological vector space method was used in the proofs. By these new 
methods we obtained an exact formula on the critical exponent /3c also in the 
case of a class of nonlinear models described by Equation ([32]). Generalization 
to higher dimensions remains a challenging open problem. 
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